Electron-ion scattering in dense mult i- component plasmas: 
application to the outer crust of an accreting neutron star 



J. Daligault and S. Gupta 
Los Alamos National Laboratory, Los Alamos, New Mexico 87545 
daligaulOlanl . gov , guptasanj ibOlanl . gov 



ABSTRACT 

The thermal conductivity of a dense multi- component plasma is critical to the modeling of 
accreting neutron stars. To this end, we perform large-scale molecular dynamics simulations to 
calculate the static structure factor of the dense multi-component plasma in the neutron star 
crust from near the photosphere-ocean boundary to the vicinity of the neutron drip point. The 
structure factors are used to validate a microscopic linear mixing rule that is valid for arbitrarily 
complex plasmas over a wide range of Coulomb couplings. The microscopic mixing rule in turn 
implies and validates the linear mixing rule for the equation of state properties and also the linear 
mixing rule for the electrical and thermal conductivities of dense multi-component plasmas. To 
make our result as useful as possible, for the specific cases of electrical and thermal conductivities, 
we provide a simple analytic fit that is valid for arbitrarily complex multi-component plasmas over 
a wide range of Coulomb couplings. We compute the thermal conductivity for a representative 
compositional profile of the outer crust of an accreting neutron star in which hundreds of nuclear 
species can be present. We utilize our results to re-examine the so-called impurity parameter 
formalism as used to characterize impure plasmas. 



Subject headings: Dense matter, nuclear reactions, nucleosynthesis, abundances, equation of state - 
neutron - X-rays; binaries, bursts 



Stars; 



1. Introduction 

The thermal conductivity of a dense multi- 
component plasma is critical to the modeling of 
the quiescent and explosive emission from the 
surface of accreting neutron stars. Much ef- 
fort has been expended in accurately calculat- 
ing the magnitude and location of crustal heat- 
ing from nuclear reactions using modern nuclear 
m odels for energetics and r e action rates such as 



( 2007 ): 



aisel fc Zdunikl dl990 l l2003h : 
Haensel fc Zdunikl (j2008f ); 



Gupta et al 



Gupta et al 



2008). However, crust thermal conductivity is 



of equal importance in determining the thermal 
structure of an accreting neutron star and the un- 
certainty in this quantity is compounded by the 
difficulties in understanding the phase diagram of 
a dense plasma with hundreds of nuclear species 
spread over several mass chains. 



For illustration of the importance of thermal 
conductivity, we now describe two phenomena 
that have attracted much attention in the observa- 
tional community: (i) superbursts and (ii) quasi- 
persistent transients. 

(i) Long-term monitoring of Low-Mass X- 
ray Binaries (LMXBs) by instruments on board 
BeppoSAX and RXTE has yielded a wealth 
of phenomena on accreting neutron stars. For 
example, superbursts, lasting roughly a thou- 
sand times as long as Type-I X-ray bursts and 
recurring on timescales of a year, have ener- 
gies ~ 10 3 times that of a typical X-ray burst. 
For the source 4U 1636-53 t hree superbursts 
were observed within 4.7 yr ( Wijnandsl 2001 



Kuulkers et al. 2004f ). These still mysterious ex 



plosions are thought to result from u nstable 12 C 
burning (jCumming fc Bildsten 2001) at column 



depths ~ 10 12 g cm" 2 (see Kuulkers et al.l (|2004h 
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for a review). The crust thermal gradient dictates 
whether a temperature conducive to unstable 12 C 
burning can exist at an ignition depth inferred 
from observations of superburst recurrence times. 
A critical ingredient in predicting the thermal gra- 
dient that a neutron star crust can support is the 
thermal conductivity of the impure crust that is 
expected when Type-I X-ray bursts drive a rapid- 
proton-capture, or "rp-process" on the neutron 
star surface. Models that treat the crust thermal 
conductivity as uncertain also have to grapple 
with the large uncertainties in core neutrino lu- 
minosity from URCA processes. In this work we 
show that if the composition of the crust can be 
accurately computed, then so can crust thermal 
conductivity. 

(ii) Another exciting application of accurate 
conductivities appears in modeling the cooling 
lightcurves of transiently accreting neutron stars. 
In the case of transients with accretion episodes 
of sufficient duration to heat the crust by deep 
crustal nuclear reactions (quasi-persistent tran- 
sients), the thermal relaxation behavior gives us a 
new observational window into the physics of mat- 



ter at v ery high density. Rec ently, IShternin et al 



(|2007t ). ICackett et ail (j2006f) have presented de- 
tailed analyses of two such quasi-persistent tran- 
sients, namely KS 1731-260 and MXB 1659-29. 
Since a critical determinant of the time evolution 
of surface flux is the thermal diffusion timescale as 
a function of depth in the neutron star, the con- 
ductivity plays an important role in modeling the 
cooling lightcurve of transients. 

To reduce the uncertainties in modeling the 
thermal profiles of persistent and transient accre- 
tors, we have begun a program of accurately com- 
puting the microphysical quantities that influence 
thermal structure. For the impure ashes from an 
X-ray burst, it is crucial to ascertain the composi- 
tion evolution of these ashes with increasing depth 
in the neutron star crust. It is upon this composi- 
tion profile that all other thermal structure deter- 
minants depend, most notably 

• nuclear heating due to non-equilibrium elec- 
tron captures, neutron reactions and pyc- 
nonuclear fusion, 

• the thermal conductivity k and 

• the neutrino cmissivity e u from Neutrino- 



Pair Bremsstrahlung (NPB). 

The compositional profile of the outer crust 
(where there is negligible free neutron abundance) 
starting with a realistic multi- component plasma 
(MCP) of XRB ashes was calculated for the first 
time using global nuclear structure and reac- 
tion rate inputs to a large reaction network in 
(jOupta et al.ll2007h . By a "composition profile" 
we refer to the changing composition from the top 
of the crust (X-ray burst ashes) all the way to 
neutron-drip. This spans four orders of magni- 
tude in density (roughly ~ 10 7 5 to ~ 10 115 g/cc) 
allowing for a large range in Coulomb coupling. 
It is this compositional profile that we employ 
in the present work to investigate the resulting 
crustal microphysics since it has tens to hundreds 
of abundant species at each depth. Thus, we have 
an ideal site to explore the physics of strongly- 
coupled multi- component plasmas (MCP). 

In this paper, we perform large-scale molecu- 
lar dynamics (MD) simulations to calculate the 
structure factor of the dense MCP in the neu- 
tron star crust from near the photosphere-ocean 
boundary to the vicinity of the neutron drip point. 
The structure factors are used to validate a mi- 
croscopic linear mixing rule (MLMR) that is valid 
for arbitrarily complex MCPs over a wide range of 
Coulomb couplings. The MLMR in turn implies 
and validates the mixing rule for the equation of 
state propertie s rece ntly proposed for MCPs by 
Potekhin etaD (|2009h . and the mixing rule for the 
electrical and therma l conductivities of MCPs by 
(Il999h . 



Potekhin et al 



These mixing rules provide 
a simple recipe to compute the physical properties 
of plasma mixtures, such as would be found in 
terrestrial experiments (e.g., inertial confinement 
fusion (ICF) plasmas.) In general, the MLMR im- 
plies a linear mixing rule for any quantity that 
is an integral of the structure factor, such as the 
neutrino pair-Brcmsstrahlung cmissivity. To make 
our result as useful as possible for the compact ob- 
ject modeling community, for the specific cases of 
electrical and thermal conductivities, we provide 
a simple analytic fit that is valid for arbitrarily 
complex MCPs over a wide range of Coulomb cou- 
plings. 

In addtion, we comment on the standard char- 
acterization of crust impurity in terms of the im- 
purity parameter Qi mp , the mean-square deviation 
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in species charge, and show that this is inadequate 
for the amorphous outer crust. We compare the 
similarities and differences charaterizing electron- 
ion scattering in the crust of neutron stars and in 
terrestrial solid and liquid metals - we also com- 
ment on the controversy regarding the behavior of 
the electron-ion scattering at the solid-liquid phase 
boundary. We use the mixing rule to compute the 
thermal conductivity of the outer crust of an ac- 
creting neutron star as a function of depth. Since 
the X-ray burst ashes (initial crust composition) 
used to determine the crust composition profile 
are obtained from a model that has accretion rate 
~ 0.1Af E dd, where M Bdd w 10 18 g s" 1 is the Ed- 
dington accretion rate for a 10 km radius neutron 
star accreting a roughly solar hydrogen fraction, 
our results for crust conductivity are immediately 
applicable to superburst progenitors. 

We restrict our study to the outer crust of the 
neutron star . In the vicinity of neutron drip, 
Gupta et al. showed that a large num- 

ber of reactions involving neutrons can remove 
compositional memory, and the complexity can 
be drastically reduced in the inner crust. We 
plan to follow up o n the nucleosynthesis results of 
Gupta et al.l ( 2008f ) with a simple model of how an 
MCP evolves in the inner crust - this will serve as 
an analytic model of crust composition evolution 
on either side of neutron-drip without invoking 
large reaction networks. Even in the scenario of 
unsuppresse d photo-dissociation as inv estigated 
recently by IShternin fc Yakovlevl (j2009h . the im- 
purity of the outer crust is preserved from X-ray 
burst ashes on the surface to the rather high depth 
of 10 11 g/cc. Between this depth and the neutron- 
drip point at ~ 4 x 10 11 g/cc, temperature-driven 
(n, 7) — (7, n) equilibria along isotopic chains 
will reduce the impurity to roughly Qimp ~ 40. 
Beyond neutron-drip (between ~ 4 x 10 11 and 
3 x 10 12 g/cc) the density-dependent electron- 
capture-dela yed neutron emissi ons, (EC,xn), as 
described in iGupta et al.l ( 2008 ). will be the pri- 
mary drivers of abundance re-arrangement and 
crust impurity will be reduced very quickly with 
increasing depth. Between - 3 x 10 12 and ~ 10 13 
g/cc pycnonuclear reactions will change the com- 
position using the few remaining nuclides to de- 
termine the compositional trajectory Finally, at 
densities greater than 10 13 g/cc proton-tunneling 
reactions will play a role, and may be instrumen- 



tal in driving the composition to a single species 



(jJones I 120051 ) without diverse metastable state 



populations spread out over different proton shell 
closures. The abrupt change in composition and 
resulting change in the phase diagram at neutron 
drip will be addressed in a later paper. Therefore 
this paper is restricted to the highly impure liq- 
uid metallic and amorphous phases of the neutron 
star crust which exist at depths shallower than 
neutron drip. 

2. Electronic properties of ionic mixtures 

The section is organized as follows: 

• First, we briefly discuss the principal chara- 
teristics of the MCPs found in the crust of 
accreting neutron star; 

• we then recall the link between the elec- 
tronic transport properties and NPB emis- 
sivity with the charge-charge structure fac- 
tor; 

• we discuss the standard treatments of mix- 
tures, namely the impurity parameter for- 
malism and the linear mixing rule; 

• finally, we compare electron scattering in ter- 
restrial metals and neutron star crusts, and 
comment on the model of electron-ion scat- 
tering at the solid-liquid phase boundary as 
proposed by lBaiko et al. ( 19981 ). 



All the formulas below are expressed in natural 
(SI) units (in particular e 2 = q 2 /Ane^ where q is 
the electronic charge.) Also, the words "nucleus 
(nuclei)" and "ion (ions)" are used interchange- 
ably. 

2.1. Characterizing the ionic mixtures 

We consider an unmagnetized, neutral multi- 
component plasma (MCP) in a volume V at tem- 
perature T consisting of a mixture of N sp atomic 
species of charge Zj and mass number Aj, j = 
1, ...,N sp . The number and mass densities of 
species j are rij and pj, respectively The total 



N S1 

and p = Y^, 7 -=i P31 respectively. The total elec- 



number density and mass density are n = ^2 



tron density is n e 



4 sp Zj n j = (Z)n, where 
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(Z) = Ylj SP x jZj is the mean ionic charge. We 
use 

Xj = rij/n (1) 

to denote the number fraction of the jth species. 
Note that the number fraction is related to the 
species "abundances" Yj as traditionally used in 
the astrophysics literature, therefore number frac- 
tion Xj differ from the usual "mass fraction" Xj. 
Thus 



n 3 = pN^Yj = pN A 



(2) 



where Aj is the mass number of species j, TVa is 
Avogadro's number and 



Xj/Aj 



Z k (X k /A k ) 



(3) 



where the index k denotes summation over all nu- 
clear species in the MCP, including free nucleons. 
The average inter-particle distance between nuclei 
is denoted by a = (3/ inn) 1 / 3 , whereas the mean 
distance between electrons is a c = (3/47T ric) 1 / 3 . 
We will occasionally use the compact notations 
Z = (Zi, . . . , Z Nsp ) and x = (x 1} . . .,s Nsp ) to de- 
note the composition. When A sp = 1 the system 
is referred to as an one-component plasma (OCP) 
and we drop the subscript referring to the parti- 
cles' species. 

We define several dimensionlcss parameters to 
characterize the system, namely: 

1. the relativistic parameter x r = p-p/m e c; 

2. the degeneracy parameter O = T/Tp; 

3. the electron coupling parameter Y c = e 2 / a c k-QT , 

where pp = hk-p is the Fermi momentum, fcp = 
{2>-K 2 n c ) 1 / d ' is the Fermi wavevector and Tp is the 
Fermi temperature of electrons. 

For the typical physical conditions of inter- 
est in this paper, the density is in the range 
10 6 g.cm~ 3 < p < 10 12 g.cm~ 3 , the temperature 
is T ~ 0.5 GK, and the number of species N sp 
varies from a few to several hundreds. As illus- 
trated in FigQ] at these densities and tempera- 
tures, x T > 1, G << 1 and a/a-Q << 1 where as 
is the Bohr radius. As a consequence, atoms are 
completely pressure-ionized and matter is a dense, 



multi- component plasma composed of bare nuclei 
neutralized by a relativistic and strongly degener- 
ate gas of delocalized electrons. This paper deals 
with the scattering of those electrons off the multi- 
component ionic mixture. 

The ions are in general non-relativistic - we 
shall also assume hereafter that the ions behave 
as classical entities, thereby neglecting the quan- 
tum nature of the ionic density fluctuations that 
may occur in the regions where T > T p = Hujp/kB, 
where uj v is the (kinetic) ion plasma frequency. 

The nuclei are "strongly coupled" in the sense 
that their properties are primarily governed by 
their mutual interactions rather than by purely 
thermal effects as in an ideal gas. 

In the OCP, electron-ion scattering is fairly well 
understood. The degree of coupling of an OCP 
is usually characterized by the so-called Coulomb 
coupling parameter Tocp = Z 2 e 2 /ak^T, the ra- 
tio of the mean interaction energy Z 2 e 2 /a to the 
thermal energy k-gT of the particles. Neglecting 
electronic screening, the OCP is known to crystal- 
lize into a bec lattice at Tocp = T c w 175. Below 
T c the OCP is in a fluid state whose properties 
vary from gas-like at small rpep < 0.1 to liquid- 
like when 50 < r OC p < T c (|Daligaurj2006i) : in the 
intermediate regime, 0.1 < Focp < 50, both po- 
tential and kinetic effects play a comparable role 
in a non-trivial way. As a consequence, models 
of neutron star crusts based on the OCP usually 
treat the electron-ion scattering as follows: 

• when r > T c , electrons travel in a bec crystal 
and electron-ion scattering is interpreted as 
the scattering of Bloch electrons off the ther- 
mal vibrations (phonons) of the ions with 
respect to the bec lattice sites. 

• when r < r c , electrons scatter off the disor- 
dered ionic background characteristic of the 
fluid state. 



Recently iBaiko et al. ( 1998 ) called into question 
this canonical picture because it yields disconti- 
nuities in the electronic transport properties at 
melting. In section l2.2.4[ we comment on this con- 
troversy regarding the behavior of the electron-ion 
scattering at the solid-liquid phase boundary of an 
OCP. 

The state of affairs for multi- component mix- 
tures is much less known, and especially for mix- 
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tures composed of N sp > 3 species as encoun- 
tered in the crust of accreting neutron stars. A 
common approach used to infer the properties of 
a mixture is to map the latter onto an effective 
OCP of coupling r o ff, thereby disregarding the 
(poorly-understood) complexities of the phase di- 
agram and dynamic properties of a true mixture. 
A suitable definition of the effective coupling is 

iV Bp 

r off = J2 x i T i 

i=i 

= {Z 5 / 3 )T C (4) 
= (Z]e 2 )/( aj k B T) (5) 



where 



is the coupling parameter for species j, a,j = 
(3/47rn J ) 1 / 3 being the average inter-particle dis- 
tance of particles of species j, and (Z 5 ^ 3 ) = 
X^ HP Xj Zj^ 3 . This choice for r c fr is motivated 
by the inequa lity (written her e for pure Coulomb 
interactions) ( Rosenfeldl 1993 ) 



10 cff " k B T ' 



(6) 



which provides a very good lower bound of the ra- 
tio of mean excess (interaction) energy per particle 
Ucx/nV to the thermal energy k B T, and therefore 
is a lower bound on the effective Coulomb coupling 
in the MCP. 

In the astrophysical literature related to dense 
plasmas, often a mixture is assumed to crystallize 
at r c ft = 175 as in the OCP. Several questions 
can be raised with regard to the onset of periodic 
ordering in a true MCP that have not been satis- 
factorily addressed to date. For instance: 

(a) Does the system crystallize in a regular, pe- 
riodic lattice as has often been assumed in 
the literature? 

(b) Do the most abundant species make a lattice 
and the less abundant species contribute to 
defects in the structure? 

(c) Is it amorphous with a liquid-like frozen-in 
structure? 

An affirmative answer to questions (a) and (b) 
above is assumed rather commonly in the liter- 
ature and consequently, beyond r e g > T c , the to- 
tal electron-ion scattering cross-section is split into 



two contributions: (1) the electron-phonon contri- 
bution due to the lattice structure, and (2) the 
electron-impurity contribution. We believe that 
the phase diagram of mixtures such as those ex- 
isting in the crust of accreting neutron stars is still 
too poorly known to make any definite statements 
such as (a) and (b) above. The local composi- 
tion may also be changed due to phase-separation 
(demixing) and sedimentation effects, which arc 
still n ot well understood a nd are being investi- 
gated (|Horowitz et al.l 120081 ). Given the complex- 
ity of the mixtures, it is, however, not unreason- 
able to think that, even when r o ff > r c , the sys- 
tem will most likely be amorphous, or at the very 
least be comprised of many defects, and that the 
existence of Bloch states is questionable. In the 
following, we shall assume that the composition at 
any depth in the outer crust of the accreting neu- 
tron star is known; in pract ice, we shall use th e 
outer crust composition from lGupta et al.l (|2007l h 
In all cases that we considered, when starting from 
an amorphous structure, the molecular dynamics 
simulations performed did not show any evidence 
for demixing or crystallization over the timescale 
of the simulations. 

2.2. Electron-ion scattering in MCPs 

We shall focus on electron-ion scattering and 
ignore the effect of ele ctron-electron (ee) sc atter- 
ing. It was shown in ( Potekhin et a.1.1 fl999h that 
ee scattering contributes to conductivity only for 
very light nuclei with Z < 6 and under non- 
degenerate conditions. For a pure 56 Fe composi- 
tion in a de generate plasma at de nsities « 10 8 5 — 
10 11 ' 5 g/cc (jPotekhin et aljfl999h showed that ee 
collisions could be neglected. 

Many of the basic formulas in t he following sec- 
tions c an be found, for instance, in lFlowers fc Itoh 
(Il976h . 



2.2.1. Thermal and Electrical Conductivities 

Electron-ion scattering impedes the transport 
of electronic momentum and energy. The result- 
ing electrical conductivity a and the thermal con- 
ductivity k can generally be expressed in terms of 
the collision frequencies v a and v K : 



m* 



(7) 
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-u f Cvt k 



-n 2 k\Tn c 
3m* 



(8) 



where m* = m e ^JT+~x* is the effective electron 
mass, vf = pp/m* is the Fermi velocity, and Cy 
is the heat capacity of the free relativistic electron 
gas. 

The quantities r a and r K describe the typi- 
cal timescale over which electron-ion collisions de- 
grade the electric and thermal currents. Their in- 
verses, the frequencies v a>K = l/r a . K , can conve- 
niently be written as the product of a "fundamen- 
tal" frequency Vq: 



4n{Z 2 )e*ni 4a 2 E F (Z 2 ) 



(9) 



multiplied by a dimcnsionlcss quantity lnA 
known as the Coulomb logarithm: 



= u Q In A c 



(10) 



In Eq.@, (Z 2 ) results form the appropriate nor- 
malization S zz (k — > oo ) — > 1 of the structure factor 
discussed below in Eq. ([T5|) . a = e 2 /he is the fine- 
structure constant and Ef = m* e c 2 is the electron 
Fermi energy. The fundamental frequency vq can 
be seen as a frequency of Coulomb collisions as- 
sociated with the transfer of momentum between 
electrons and ions; it can be written as 



vq = n 



Air 



(Z 2 )ie 2 



VF, 



111) 



in which the bracketed term can be interpreted 
as the cross-section for momentum transfer in a 
binary collision between an electron of velocity 
and charge — e and a massive ion of charge 
(Z 2 )?e. The dimensionless Coulomb logarithm 
contains all the fine details of electron-ion colli- 
sions in a plasma, such as charge screening, col- 
lective modes and quantum diffraction effects. In 
general, lnA CT ^ lnA K since the two transport co- 
efficients measure different mechanisms, namely 
momentum and energy transport. However, when 
the dominant source of electron-ion scattering is 
clastic, as in the disordered fluid phase or the high- 
temperature crystalline phase, the two Coulomb 
logarithms are equal and the Wiedemann-Franz 
relation 



k/ct 



(ir 2 k 2 B /3e 2 )T 



(12) 



is satisfied. 

A calculation of the Coulomb logarithm that 
self-consistently treats all the many-body effects 
occuring in a multi-component plasma is in gen- 
eral rather complicated. However, because of the 
high electron density in neutron stars, the typical 
electron kinetic energy is very large compared to 
the electron-ion interaction energy and the average 
electron-ion correlations are weak. Accordingly, a 
good approximation can be made by considering 
the ions and electrons in the plasma as two weakly 
interacting subsystems, namely an MCP neutral- 
ized by a homogeneous negatively-charged back- 
ground and an homogeneous interacting electron 
gas (i.e. relativistic quantum electronic jellium.) 
Electrons scatter off the electrostatic potential cre- 
ated by the ionic charge distribution, 



n 2 (r,t)=^Z J ^<5(r-R i (t)), 



(13) 



i=i 



z=i 



where Rj(f) is the position of the l-th. nuclei of 
the j-th species at time t, and Nj = njV is the 
number of particles of species j. In the (first) Born 
approximation^ of the electron-ion interaction and 
using the fact that the el ectrons are much heav ier 
than the ions, we obtain ([Flowers fc Itohlll976t ): 



lnA„ 



2A T 



dkt 



v(k) 



e e (k) 



k 
2k~F~ 



(14) 



with v(k) = 1/k 2 and e e (k) is the static dielectric 
function of the electron jellium. 



1 Extensions to higher-order approxim ations are difficult. An 
approximate treatment proposed bv lYakovlev et at I <1987fl 
prescribes usage of the exact binary cross-section instead 
of the Born approximation in the integrand of Eq. lldl . 
This introduces an additional factor in Eq, l|14j l. namely 
R(k) = o-(fc)/<7g(fc) where a is the exact differential cross- 
section for a momentum transfer k and ag is its Born ap- 
proximation. A similar procedure can be applied verbatim 
to our calculations since the MCP structure factor is not af- 
fected. Re c ently, a similar averaging procedure was used by 

lltoh et al. 

I j2008t) , however the title of that paper is some- 
what misleading since it is not actually the second-Born ap- 
proximation that is performed therein. Such a calculation 
would involve the higher-order correlation functions of the 
ionic density Eg . II 13t . However, these higher-or der correla- 
tions are noticeably absent in the equations of lltoh et al.l 
(2008) for the transport coefficients. 
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Equation (fT4| can be interpreted as follows. 
The Coulomb logarithm lnA ajK sums all the con- 
tributions of collisions characterized by momen- 
tum transfer k. The range of integration is limited 
to 2fcp because the electrons are degenerate and 
only those electrons near the Fermi surface take 
part in exchanges of momentum and energy. The 
electrons behave as if they were scattered indepen- 
dently of each other by the electronically screened 
electron-ion potential. Because the electrons are 
much lighter than the nuclei, the energy exchanges 
fko in inelastic collisions occur in a range over 
which e e (k,u>) « e e (k, 0) = e e (k). For relativis- 
tic, degenerate electro ns, an analytica l expression 
for e e was obtained bv ljancovicil (|l962f ) within the 
Random-Phase Approximation (see Appendix [Bj . 

The bracketed term, which describes the kine- 
matic suppression of backward scattering of rel- 

tween unity for x T — > and [l — 
x T —> ac. 

The many-body physics of the ionic subsystem 
is encapsulated in the function Sa- tK (k). S a ^(k) is 
a functional of the spectrum of charged fluctua- 



ativistic electrons, 



1 



varies be- 
(k/2k F ) 2 ] for 



tions 5n z (k, t) 
system 

S zz {k,u>) = 



n z (k,t) — (n z (k)) eq in the ionic 



(15) 



iV(Z 2 



dio 
2^ 



(5n z (k,t) 5n z 



-k,0)) eq e* 



where ( ) cq denotes the average of a canonical en- 
semble at temperature T. In other words, S a ^ K (k) 
self-consistently takes into account all of the sep- 
arate electron-ion relaxation mechanisms at the 
level of the Born approximation; it is really an 
umbrella term for electron-ion, electron-impurity, 
electron-phonon scattering and includes multi- 
phonon processes. We have, from lFlowers fc Itoh 
i|l976h . 



S a {k) = / S zz (k, Lu)zn(z) z = hw/kBT 

J — oo 

S K (k) = J S zz {k^)zn{z) l+l^-- 



where n{z) = 1/1 — e z . 

In the special case of a periodic, crystalline 
structure, the term (n z (k)) oq peaks at values of 



the wavevector k commensurate with the period- 
icity of the lattice and substracts out the elas- 
tic Bragg peaks: electron-ion scattering involves 
Bloch electrons and Eq. {HP is obtained by approx- 
imating their behavior by means of free particle 
wave functions. 

When the scattering is predominantly elastic in 
nature, as is expected in the fluid and solid phases 
at T > T p or in a disordered (non-periodic) lattice, 
we have 



S zz (k,u)) 



(16) 



OO 

1 



nV(Z 2 ) 
S zz (k) 



(6n z (k,0) 5n z (-k, 



(18) 



where S zz {k) is the (charge-charge) structure fac- 
tor. This leads to 



lnA CT , K = 

2fc F 



dkk : 



v{k) 



1 - 



1 



k 
2k^ 



(19) 
S zz {k\ 



In this limit, a and k are insensitive to the details 
of the spectrum of ionic charge fluctuations. These 
details (ion collective modes and single-particle ef- 
fects) are present implicitly but in an integrated 
form via the sum- rule (|18(1 . 

Until now, we have assumed that the ions are 
point particles. In this "point-charge" approxima- 
tion, the only species-dependent term in the in- 
tegrand of Eq. (fT4j) is the charge-charge structure 
factor S zz (k). In order to allow for finite nuclear 
size, one can replace the density (|13p with 

JV Bp N 3 

n z (k,t)=Y,Z ] F ] (k)J2e- lk - Rlit) - (20) 
j=i i=i 

where Fj(k) is the nuclear form factor, which re- 
flects the charge distribution within a nucleus of 
species j. For instance, assuming uniform proton 
density in the nuclei, 



Fj(k) = (smfrRj) - kRj cosCfcJ?,)) (21) 

where Rj the charge radius of the nuclear species j 



which has mass number A , ; we use Rj = 1.15 A 



1/3 



— „„„„„ " t y 3 

fm as in iKaminker et al (1999). Direct MD cal- 
culation of the structure factor S zz (k) with and 
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without the nuclear form factor (|21j) shows that 
the finite nuclear size corrections are completely 
negligible in the outer crust of interest in this pa- 
per. (Note that for instance (F^k)) 2 = 0.9766 
at the peak of S zz (k), i.e. at (k/2kp) ps 0.33, 
for a typical heavy nucleus having mass number 
Aj ps 64 - 100 or A 1 / 3 ps 4 - 4.6). Accordingly we 
neglect finite nuclear size effects in the remainder 
of this paper. 

2.2.2. Neutrino-pair Bremsstrahlung (NPB) emis- 
sivity 

In NPB, the electrons are accelerated by 
the Coulomb field of the crust as in photon 
Bremsstrahlung but with the emission of a neutrino- 
antineutrino pair instead of a photon. In the first 
Born approximation, the neutrino emissivity is 



be used to compute both k and e„ for an MCP of 
arbitrary complexity using the microscopic mixing 
rule (MLMR) described in section[l]which involves 
only the OCP structure factors and species abun- 
dances. For higher values of T > 225 in the neu- 
tron star crust we present a simple model for the 
behavior of the Coulomb logarithm that reduces 
the evaluation of thermal and electrical conduc- 
tivity to an analytic formula incorporating crust 
impurity in a very simple way. We do not present 
calculations of the NPB neutrino emissivity in this 
paper, rather we point it out to demonstrate that 
all such "integral" determinants (i.e. summations 
over momenta) of crust thermal profile are related 
to the crust composition via the quantity S zz {k). 

2.2.3. The standard treatment of mixtures in the 
crust 



8ttG f (Z 2 



567** (fcBT)6lnAt 
with the Coulomb logarithm for NPB, 



In A f = 



where k 



he 

for 

2/,-f 



2 /, T 



dkk c 



dk t fcjf / dk T 



v(k) 



e(k) 



o 



v(k) 



e(k) 



a(k)S zz (k) 



k t 



r , Gf is the Fermi weak cou- 
pling constant, C\ ps 1.675 from the vector and 
axial- vector constants of the weak interaction and 



^22) F° r simplicity, we drop the a, k and v subscripts 

in what follows. 

The most common approach of electron-ion 
scattering in the neutron star crust relies on the 
' S zz (k)B,T(kt, kpp^ e ^ that the MCP is fluid in the upper part 
of the crust where r e ff < 175 and consist of a 
solid, i.e. a crystalline lattice with impurities, 
the lower part where T e g > 175. Accord- 
ingly, in the solid phase, the collision frequency 



a(k) = 



1 - 



2 k- 



fr- 



ill 



(2^) 



The' function i?r 



describes the thermal broadening of the electron 



v = z^e-ph + ^c-imp is split into the sum of the 
electron-phonon and the electron-impurity scat- 
tering contributions; in the liquid state v = is C i 
where v e \ is the electron-ion scattering inverse re- 
laxation time discussed in the next section. 



Fermi surface and from Eq.(20) of lHaensel et al 



( 19961 ) this can be approximated as unity for T < 
0.6 GK typical of models of the neutron star crust 
at accret ion rates characteristi c of supcrburst pro- 
genitors ( Kaminker et al 19991 ). 

Thus, in both Ea.([Ii]l and Ea.([2"3"|) the only 
dependence on the nuclear species present (i.e., 
on the composition) enters via S zz {k). Therefore, 
to tabulate fits to both k and e v for every nu- 
clear species prior to application of a mixing rule 
i s extremely cumberso me (note that the papers 
(|ltoh et al.ll2008lll996h tabulate fit coefficients for 
only 11 species for conductivities and NPB emis- 
sivities respectively) and also unnecessary from a 
computationa l standpoint. The fi ts to OCP struc- 
ture factors in I Young et al.l <| 1991! ) are sufficient to 
cover every one of the N sp > 100 species encoun- 
tered in the neutron star crust. They can also 



The splitting v 



tially introduced in Flowers fe Itoh ( 19761 ) to deal 
with periodic crystals with a small fraction of 
the lattice sites occupied by impurities. The 
idea is illustrated in Fig.(j2]). The case of bi- 
n ary mixtures (Asp = 2) was discussed in detail 



Itoh fc Kohvamal (|l993T) for a small concentra- 
tion of impurity ions (say species "2"): the con- 
tribution of the impurity to electron-ion scatter- 
ing is equal to its concentration times the binary 
electron-ion collision frequency between an elec- 
tron and a "residual" ion of charge (Z2 — Z\). 
This idea was extended to deal with more com- 
plex MCPs containing A sp > 3 species, in which 
case the system is treated as a perfect crystal of 
charge (Z) plus the contribution of uncorrelated 
impurities of charge (Z\ — (Z)). 



S 



Within this approximation and defining the im- 
purity parameter as 



Qimp ~ ^ Xj{Zj — (Z)) 2 



(24) 



the relaxation frequency can be written as 

^imp = ^O.imp l n A; mp (25) 



where 



and 



4a 2 e F Q imp 
3irh (Z) 



(26) 



In A, 



imp 



2fe P 



dk¥ 



v{k) 



i 



1 + x 2 \2k 



lnA; mp corresponds to Eq. (f2"0")) with S zz set to 
unity to correspond to independent (uncorrelated) 
electron-impurity collisions. 

While the impurity parameter formalism was 
developed to treat small impurity concentration 
in an otherwise almost perfect crystal lattice (see 
Fig. ([2]), it has also used to characterize electron- 
ion scattering in complex fluid and amorphous 
MCPs when Qi m p is large, say Q 



imp 



100. Such 

is usually the case when electron-impurity scatter- 
ing is used to set a lower bound on conductivity 
for the entire crust. However, since A; mp is of the 
order of magnitude as A, we have 



v/v- 



In A 



imp 



Q'\m\> In A-imp 



jimp 
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(28) 
(29) 



over the whole outer crust as shown in FigfTJ 
the impurity parameter can lead to conductivi- 
ties much higher than the those that would be 
obtained using the actual structure factors cm- 
ployed in Eq.(|20|). We thus anticipate that Qi mp 
is inappropriate for the outer crust (impure liq- 
uid metal approaching heterogeneous amorphous 
solid at high density), and instead we shall justify 
the usage of (Z 2 ) from its presence in our expres- 
sion for the fundamental frequency vq in equation 
Eq.©. The differences in crust conductivity that 
would result from these differing prescriptions are 
discussed in section [5] 



(27) 



In the liquid state, Potekhin et al. ( 19991 ) pro- 
posed a different approach to deal with mixtures. 
In this formalism one writes, 



In A 



1 



x- J Z, 2 lnA OCP (r,) (30) 



where I\ 



z 5/3 r„ 



Xj refers to the number frac- 



■ J ~ " j 6! J 

tion as defined in Eq.fl]), and lnA OCP (r) is the 
Coulomb logarithm (|2T)|) for an OCP at coupling 
r. This prescription for mixtures was also rec- 



omme nded recently in the w orks of iCassisi et al 



(|2007t ) and lltoh et all (|2008l) and we corroborate 
its efficacy in section |4j This presciption (|30l) . 
when substituted in Eq. (TlT)]) . leads to 



v = Vq- 



[Z 2 lnA OCP ) 



(31) 



In sections @] and EJ we propose and validate 
with MD simulations a mixing rule for \lnA that is 
slightly modified from Eq. (|3"0)) . We also show that 
the impurity parameter formalism is inappropriate 
at large Q imp . 

2.2.4- Discontinuities in transport coefficients at 
crystallization 

The previous expressions tell us that the con- 
ductivities and NPB emissivity depend on the 
scattering of electrons off the fluctuations 8n z (k, t) 
in the ionic charge density about the average value 
(n z (k)) eq . However, this average charge density is 
dependent upon thermodynamic phase of the sys- 
tem: 

• If the system is fluid, the averaged charge 
density (n 2 (k)) oq = n e (5(k) is uniform; 

• If the system is in a crystalline solid state, 
i.e. comprised of a periodic (Bravais) lat- 
tice with the same physical unit, elec- 
trons arc in Bloch states and (n z (k)) = 
e -2W(k) n ^y § k K depends on k, where, 
assuming A^p = 1, K are reciprocal lattice 
vectors and W(k) is the Debye- Waller fac- 
tor: Bloch electrons do not scatter off the 
static, wavevector-dependent charge fluctu- 
ations (Bragg peak). 

Because of this fundamental difference, the for- 
mulae of the previous section predict significant 
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(by a factor of 3-5) discontinuities of the electric 
and thermal con ductivities a t the melting point 
(jltoh et al.l fl993h . Recently. Iriaiko et all (|l998h 
argued that this jump in the OCP electronic trans- 
port properties upon crystallization is not physi- 
cal and that the electron transport properties in 
strongly coupled plasmas should be fairly insen- 
sitive to the thermodynamic state of the plasma. 
Accordingly, they advance that strongly coupled 
Coulomb liquids are characterized locally by the 
same long-range order as present in the crys- 
talline solid and although the long-range order 
does not persist forever as in a solid (the liquid 
state is fluid), electrons keep traveling in those 
local crystal-like structures without any signifi- 
cant degradation of their mean velocity. Thus, 
Baiko et al.l (|1998l ) advocate that, in the liquid 
state, the contribution corresponding to elastic 
scattering off the incipient crystalline structure 
must be substracted from the static structure fac- 
tor. As intended, this prescription removes the 
jump in the conductivities at the solid-liquid tran- 
sition - this prescription was later used in a se- 
ries of papers to redefine practical fitting formu- 
las for the el e ctron conductivities and NPB (e.g, 
Cassisi et all (120071) and references therein ; also 
Horowitz et al.l (|2008l )V 



The prescription of iBaiko et al.l (|1998f ) contra- 
dicts our current understanding of solid and liquid 
metals. We instead agree with Itoh et al. ( 2008h 
since we believe that the disordered positions and 
the oscillating and diffusive movement of the ions 
cannot be ignored in a liquid. 



The arguments and conclusions of I Baiko et al 



(|1998j ) contradict both experimental findings and 
the contemporary theoretical understanding of liq- 
uids: for example, it has been known for a long 
time that the electrical resistance of most metals in 
the liquid state just above their melting points are 
abo ut 1.5-2.3 times larger th an those of solid met- 
als ( Cusak fc Enderbv 1960l ). Even the simplest 
metals, whose electronic structure do not change 
appreciably upon melting, show a discontinuity in 
transport properties at melting, e.g. for sodium 
a s /a\ ~ 1.45 and for aluminum <j s /u\ ~ 2.2. In- 
terestingly, similar to Baiko et al.'s approach but 
with the opposite purpose in mind, early theoreti- 
cal attempts to explain the increase of electrical re- 
sistivity upon melting relied heavily on ideas that 
had proved successful for solid metals. These ap- 



proaches overlooked the effect of disorder inherent 
to liquids: for instance, in the "quasi-crystalline" 
models, the local coordination just above the melt- 
ing point was treated very similar to that which 
prevails in the solid phase just below. It was also 
believed that, as in a solid metal, the resistivity of 
liquid metals could be divided into two parts (a) 
a thermal term, proportional to the mean square 
amplitude of vibration of the ions, and hence de- 
pendent on the temperature T, and (b) a residual 
term, independent of T due to deffec t s (e.g ., va- 
cancies, dislocations). In 1934, Mottl ( 19341 ) pre- 
sented an elegant theory to reconcile these solid- 
like picture to the increase of clcctical resistivity 
on melting in terms of the entropy change on melt- 
ing per particle As = (S\ — S s )/nV, 



As = 1.5&B In (cr s /(7i) 



(32) 



Using the latent heat of fusion to evaluate As, 
Mott's model yields surprisingly good values of 
<j s /(j\ (1.68 for Sodium instead of 1.45.) It was 
later realized that Mott's assumptions do not hold 
in the light of other facts. Hence Mott's theory 
assumes that all the melting entropy is connected 
with the thermal agitation of the ions and ignored 
the configurational entropy arising from disorder 
in the rest positions about which ions are instan- 
taneously vibrating (the so-called "cage effect".) 
Unfortunately, this argument does not hold: it is 
indeed known that for all simple metals, the exper- 
imental As is approximately a universal constant, 
As « 0.8 ± 0.1 (rms), and is chi efly due to th e 
change in configurational entropy (|Wallacdll997l ). 
The entropy change due to the thermal agitation, 
the anharmonicitics in ionic motions and the elec- 
tronic structure are responsible for the small scat- 
tering in the universal value. As illustrated in ta- 
ble [TJ similar findings apply the OCP over a wide 
range of electronic screening, namely As ~ 0.85. 

Finally, we would like to add another argument 
in disfavour of Baiko et al.'s prescription. If crys- 
tallites do indeed exist in the liquid phase as as- 
sumed by IBaiko et al. (1998). they must be ex- 
tremely small, smaller than required for a local 
band structure to be established, because the pres- 
ence of nuclei would make it impossible for a liquid 
to supercool. However, all liquid metals can be 
persuaded to supercool through as much as 20% 
of th eir melting temperature (jTurnbull fc Cech 



1950t ) and MD simulations of the OCP show the 
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same trend (lDaligaultll2007t ). 

Interestingly, for the multi-component systems 
found in the crusts of accreting neutron stars, the 
jump may be much le ss important than suggested 
by the calculations of lltoh et al.l (|1993l ). However 
this is likely not due to the reason proposed by 
Baiko et al. but because of the lack of long-range 
order for those complex mixtures and the high de- 
gree of impurity-like disorder that those systems 
possess. 

2.2.5. Terrestrial liquid metal compared to the 
neutron star crust plasma 

In this section we continue the comparison of 
electron-ion scattering as occurring in terrestrial 
liquid metals as opposed to the crust of a neutron 
star. To make this comparison as quantitative as 
possible, we shall discuss the different components 
of the Coulomb logarithm in Eq. {T3]| . While this 
expression applies to both terrestrial liquid met- 
als and to the crusts of accreting neutron stars 
(under the assumption of the first Born approx- 
imation), there are several interesting differences 
between both systems which we enumerate below: 

(a) While in the former the number of species is 
typically of order unity (N sp = 1 in simple 
metals; N sp = 2 — 3 in alloys), the num- 
ber of species in the latter can be several 
hundred (see section [33]) Nevertheless, as 
shall see, the charge-charge structure factors 
S zz (k) is similar in both cases, i.e. char- 
acterized by damped oscillations originating 
from the short-range order of strongly cou- 
pled systems. The heights of the peaks in- 
crease with the strength of the coupling pa- 
rameter r e gf and therefore with depth. 

(b) While in terrestrial liquid metals electrons 
are non-relativistic, electrons in an accret- 
ing neutron star crust are relativistic. The 
main effects of relativity are encapsulated 



1 - 



l+x- 



(Jl-X 

\2k F J 



which varies from 1 



when 



and 1 - (k/2k F ) 2 when 



oo and describes kinematic suppresion of 
backward scattering of relativistic electrons. 
Relavity also modifies the electronic screen- 
ing of the ions (see Eq. (|B2[) . The over- 
all effects of relativistic corrections is il- 



lustrated in Fig.© and in Fig.lgJ). Fig- 
ure ([3]), shows the integrand of Eq. (fT4|) for 
both relativistic and non-relativistic elec- 
trons (x r = in Eq.(Q3])). Relativistic 
effects are noticeable in large momentum 
transfer collisions (close encounters) and, as 
a consequence of the suppression of back- 
scattering, result in lower scattering cross- 
sections and consequently higher conductiv- 
ities. Fig. ([D) compares the relativistic vs 
nonrclativistic Coulomb logarithm as a func- 
tio n of depth us i ng the compositional profile 
of iGupta et al.l (|2007T ) discussed in section 

(c) While in the neutron star crust, ions arc 
fully stripped and all electrons are delocal- 
ized electrons and participate in the elec- 
tronic conduction, n e = (Z)n with (Z) ~ 40, 
most electrons in a liquid metal are bound 
to the nuclei, n e — (Z c g)n where Z e g is the 
effective charge of the ions and Z c s = 1 — 3. 
As a consequence, in Eq. (fT4"|) , Vi e (k) must 
in principle be replaced by a pseudopoten- 
tial that mimics the effect of bound electrons 
on free-electron ion scattering (this is a ma- 
jor difficulty in condensed matter physics). 
In contrast, for neutron star crusts, Vi e is 
known: it is the pure Coulomb potential, 
slightly corrected with a form factor to de- 
scribe the finite size of the nuclei and very 
close electron-ion collisions. 

(d) A final effect of major importance arises 
from the magnitude of the electron Fermi 
wavevector fcp = (37r 2 n e ) 1//3 , which increases 
like Z^ff in terrestrial liquid metals as op- 
posed to (Z) 1 / 3 in the neutron star crust. 
Therefore the ranges of integration < k < 
2k-p in Eq. (fT4f greatly vary between the two 
cases. The situation is illustrated in Fig.©: 
the limit of integration 2fcp lies just to the 
left of the main peak in S(k) for mono- 
valent liquid metals, just to its right for 
divalent ones, and progressively further to 
the right for metals which have a valency 
greater that two and three. In neutron stars, 
2fcp can be quite large (and increases with 
depth), i.e. large momentum transfer colli- 
sions (close encounters) become much more 
important in determining the electron-ion 
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scattering cross-section. 

3. Molecular Dynamics (MD) Simulations 
of the neutron star crust 

In this section, we discuss the MD simulations 
we performed in order to calculate the structure 
factors Szz(k) of MCPs representative of neutron 
star crust and to validate the linear mixing rule 
discussed in section |U 

3.1. Input : a representative crust compo- 
sition profile 

To base our calculations of crustal k and e v on a 
realistic pre-neutron-drip crust comp ositional pro 



file we have used the compositions of lGupta et al 



(|2007l ) . Our choice was determined by a number of 
factors. First, the diversity of species in the start 



i ng co mposition (XRB ashes from ISchatz et al 



(|200ll ) is characterized by the "impurity param- 



imp 



100 



an 



eter" which has a high value Q 
ideal test-bed for validation of a mixing rule in 
an MCP. In the outer crust, prior to the onset 
of neutron reactions, electron captures preserve 
the high impurity. Only beyond neutron-drip 
in the inner crust, as demonstrated in the SEC- 



nucleosynthesis process described in iGupta et al 



(|2008f ), do electron-capture-delayed- neutron emis- 
sions rearrange abundances between very diverse 
mass chains. Finally, the accretion rate closely 
match those expected of X-ray superburst progen- 
itors (~ O.lMEdd ~ 10 17 g s -1 ) and therefore we 
intend the numbers we compute to have practical 
relevance to the neutron star crust modeling and 
X-ray observational communities. 

The mixtures we have simulated discard the less 
abundant species which would have had less than 
10 particles in the simulations. Details for each of 
the compostions used at depth are given in table 
[2] and the spread in Z is illustrated in Fig.© for 
two compositions. 

We now use the composition profiles at various 
neutron star crust depths to validate the linear 
mixing rule. We adopt a fiducial temperature T = 
0.5 GK, and vary the density as in Table [2] Thus, 
changes in r e g with changing crust depth arise 
only from composition evolution and a decrease in 
inter-ionic spacing . At each depth a large number 
of mass chains A ss 20 — 100 enter the composition 
as expected from X-ray burst ashes evolving in the 



outer crust at depths shallower than the neutron- 
drip point. 

3.2. The MD Procedure 

We employ molecular dynamics (MD) simula- 
tions to calculate the structure factor (fT8|) at vari- 
ous depths of a neutron star crust. In order to deal 
with the impure multi-component ionic mixtures 
typical of neutron star crusts, we have developed a 
code that allows for the simulation of large charged 
systems (tens to a few hundred nuclear species) 
over the long time scales required for the equili- 
bration of such impure MCPs. The code is based 
on a parallelized implementation of the Particle- 
Particle-Particle-Mesh (P 3 M) algorithm with pe- 
riodic boundary conditions; high resolution for in- 
dividual encounters is combined with rapid, mesh- 
based, long range force calculations. It simulates 
the classical dynamics of a mixture of nuclei de- 
scribed as point particles interacting through a 
screened Coulomb interaction. Thus the interac- 
tion potential between two nuclei of charge Zi and 
Zj is modeled by the screened Coulomb (Yukawa) 
potential, 



Vij (r) 



ZjZje k sc r 



(33) 



where fc sc is the inverse screening length. The 
P 3 M algorithm can handle any value fc sc > 
to very high accuracy. In practice, we use the 
inverse of the relativistic Thomas-Fermi inverse 
screening length fc r TF for fc sc given by k r Tpa e = 
0.185(1 + X- 2 ) 1 / 4 . 

Several technical details of the MD calcula- 
tions (simulation length, number of particles, time 
steps, etc.) are collected in appendix [A] Here 
we briefly outline the procedure used to calculate 
S zz (k). Given a composition at a certain depth, 
the particles are initially randomly distributed in 
a cubic box of volume V — L 3 with a Maxwellian 
velocity distribution at the prescribed tempera- 
ture. After a long equilibration phase, a simula- 
tion is performed in which, at each time step n, the 
charge density nz (k; n) as in Eq. (fT"3]) is calculated 
for values of k commensurate with the simulation 
box size L, namely 



Ikll < 2fa 



(34) 



where n x y z are natural integers. The charge- 
charge structure factor (fT5)) is calculated at the 
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end of the run using 

N. 



Validation of the Linear Mixing rule 
and Practical Fit 



S zz (k) 



N, 



— Yl -^K(k,n)| 2 . (35) 4.1. Validation of the LMR 



n=\ k , | |k| |=fc 



3.3. Molecular dynamics structure factors 

MD results for the structure factor S zz (k) 
and corresponding to the crustal compositions of 
Gupta et all (|2007h listed in table ©, all at tem- 
perature T = 0.5 GK, are shown in Figs. (05]) over 
the range < k < 2kp. Additional data obtained 
from the MD simulations to test lower tempera- 
tures near the photosphere-ocean boundary and 
potentially high temperatures at different depths 
are provided in the Appendix [Cl 

The shape of the charge-charge structure fac- 
tor S zz (k) is typical of the structure factor of liq- 
uids. The first peak, which systematically oc- 
curs at k/2kp ~ 0.33 and which has a height 
increasing with the effective Coulomb coupling 
r o ff, reflects the existence of a dominant short- 
range order of the particles in real space. This 
is illustrated in Fig.© that shows different pair- 
distribution functions g a b(f) for several pairs of 
species a and b. The sharp decrease in the pair- 
distribution functions at small separation, origi- 
nating from the repulsion between like charges, is 
responsible for the subsequent maxima and min- 
ima of 5 zz (fc), whose oscillation is strongly damped 
as k increases. Eventually, at large k, S zz {k) ap- 
proaches unity, in accordance with the normaliza- 
tion in Eq (fT5j) . At small k, S zz (k) probes the long- 
wavelength static fluctuations in the charge den- 
sity and approaches zero (as k 2 ) in the k — ► limit 
as a consequence of charge neutrality and perfect 
screening. 

All the simulations at density larger than ~ 
10 10 g.cm -3 (compositions #9-17 of table [2]) and 
characterized by r e ff > 260 show pronounced 
short-range structure as seen in the shoulder of the 
second peak of S zz {k). Nevertheless, within the 
timescale of the simulations, we did not witness 
any clear-cut phase separation or crystallization. 



A common approximation for the excess (non- 
ideal) energy U cx and free energy F cx of strongly 
coupled ion mixtures is the (empirical) Linear 
Mixing Rule (LMR). In terms of the energy per 
particle and per unit of ksT, u CK = Uc^/NksT 
and / ex = F cx /Nk B T, the LMR claims that 



Wex(r e , 1?) 



/ex(r e ,xO 



LMR 

N ST3 



(r e ,-*) 



,OCP 



'3 "ox 



(Ii) 
3=1 

/c L x MR (r e ,^) 

x j/ex (IjO 

3=1 



(36) 



(37) 



where T,- 



- Z^ 

excess free energy of an OCP at coupling Tj and 



ZV T e , and m™ r and / e L x MR are the 



=N S , 



"x* = {xi} i= ^ sp is the composition vector of species 
number fractions (see Eq.([T])) for the MCP. 

To the best of our knowledge, the mixing rules 
(f36|) and (f37|) were obtained empirically and there 
is no rigorous "derivation" of them. Equation (|37|) 
was shown to be very accurate for binary mixtures 
with both rigid (jPeWitt fc Slatervll2003l) and po- 
lariza ble backgro und electrons (|Chabrier fc Ashcroftl 
1990l ). Recently. Fotekhin et all (|2009h proposed 
to extend the LMR to calculation of the equation 
of state (EOS) of multi-component mixtures. The 
results derived in the following valida te this recent 
prescription of IPotekhin et al" ( 2009i ). 

The LMRs (j36)) and ((37)) suggest a more funda- 
mental, microscopic mixing rule that directly in- 
volves the structure factors. To this end, we recall 
that the excess energy of a screened MCP (i.e. 
with ions interaction via the screened Coulomb 
potential ([55)) ) can be obtained from t he charge- 
charg e structure factor S zz (k) using (jRosenfeld 
19931 ) 



w ox (r e , it) = 
r e (z 2 ) 



dk 



k 2 + k 2 



■[S zz (k)-1] 



(38) 
(Z 2 ) 



(ZV/3 



The exact equation (|55|) together with the accu- 
rate approximation Eq. ([5T]) suggest the following 
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Microscopic Linear Mixing Rule (MLMR) approx- 
imation for the MCP structure factor 5 zz (/c), 



<7 (h-V ~ QMLMR 

'-'zzl^n-) 1 e, X- ) ~ "-> zz 

(Z) 1/3 



(fc;r e ,3?) (39) 



where Sp cp (fc;r) is the structure factor of the 
OCP comprised of species j at coupling T. In- 
deed, introducing Eq. (f3"9")) in Eg. ([55]) , we obtain 
the LMR (Ml) and (J3ZJ). By introducing Eq.® 
in Ea. (fl"4)) . we obtain the LMR for the Coulomb 
logarithm, 



In A 



(z 2 ) 



E 



tjZ) ' In A ( " '': 1',. Z,. 



where In A OCP (rj, Zj, x r ) is the Coulomb loga- 
rithm for an OCP of charge Zj at the coupling 
Tj and relativistic parameter x T (the explicit de- 
pendance on the charge Zj comes from the upper 
limit 2k F = 2(3ttV) 1 / 3 = 2(3ir 2 Z ] n ] ) 1 / 3 of the 
integral in Eq.fTip). 

The validity of the MLMR is confirmed by 
our MD simulations. The structure factors for the 
compositions (1) to (8) of tabled corresponding to 
8.6 • 10 6 < p < 4.8 • 10 9 g.cm- 3 and 70 < L < 140, 
are shown in Fig. (JT|) and Fig.@. The blue lines 
show the results obtained using MD simulations 
and the red lines show the results obtained with 
the mixing rule Eq. (j39|) . The latter was evaluated 
with the structure factors for the one-component 
plasmas obtained b y solving the so-cal led HNC 
equations (see e.g. I Young et al. ( 199l[ )) For all 
the compositions, Figs.(p} and ([2]) show very good 
agreement between the MD and the mixing rule 
values. In each case, the position of the peaks is 
well reproduced and the mixing rule overestimates 
their heights by at most 10%. (The MLMR is also 
illustrated in the figures of appendix [C]) 

When used in Eq. (J38J) or Eq. (fT4"]) . however, 
these differences between the MD and MLMR 
structure factors barely affect the value of the 
Coulomb logarithm, and in turn the electronic 
transport properties are unaffected. This is illus- 
trated in Fig. |4] that shows the sum lnA(A") = 

So dkk3 38 2 f 1 ~ r& { afc) 1 with < 



rule results. Thus the MLMR provides a very 
accurate and convenient way of calculating the 
electronic transport coefficients for complex ionic 
mixtures in the liquid phase. More generally, all 
quantities that are functionals of the structure 
factor such as (f38|) . the electrical and thermal 
conductivity and even the NPB emissivities can 
be obtained to high accuracy from the MLMR. 

4.2. Practical fit to the Coulomb loga- 
rithm of a mixture 

Having validated the LMR, we now provide a 
simple practical fit to the Coulomb logarithm of 
an arbitrarily complex liquid or amorphous multi- 
component plasmas. The fit is valid for any ef- 
fective Coulomb coupling Tj > 10 and relativis- 
tic parameter range 1 < x T < 1000. For the 
neutron star crust applications in the amorphous 
state, once the composition as function of depth is 
known, the conductivity can be calculated very ac- 
curately (when compared to numerical MD simu- 
lations) without resort to complicated tabulations 
of fit coefficients or decomposition of conductivity 
from different scattering processes. 

Our prescription is as follows: 

ln A - £ ^ Z T ^ A ° cp ( r - ( 41 ) 

where Tj = Z 2 e 2 /a k B T and lnA OCP is fitted by 
the expression 



\nA OCP (T^Z^Xr) 
~ C(Zj, Tj, x r 



l + o.mrj + o.ooooir^ 2 

l + aix^Tj+bix^T 2 



ln A imp (- 



In Eq.jZZP, 

ln A imp (a; r ) 
1 



1 + 396a; r 
0.466 + 394x r 



'1 + 4/3V ) ln (1 + q- 2 ) -f3 2 - 



(43) 

1 + /3V 
l + q 2 



K < 2fcp obtained with the MD and the mixing 



is obtained by setting e(k) — 1 + k 2 TF /k 2 ii 



Vf x t 
fc rTF 0.048196 



Ea.(|27j). where 
Q 



2fo 



VP 
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fc r TF is the relativistic Thomas-Fermi wavevector, 
and the functions a(x r ) and b(x T ) are given by: 

1+a ^ 1/3 z r <5 

(44) 



a(x r ) 



a2+a-3X T 



2.22 x r > 5 



This fit was obtained to match the OCP struc- 



ture factors of lYoung et al.l (|199ll ) for the range 



5 < Tj < 225. However, assuming that the system 
stays amorphous (no periodic lattice structure), 
the fit has not restriction on Tj as demonstrated 
in FigEE! 



b(x r ) 



X r < 10 



(45) 



0.087 x T > 10 



with the parameters a n and b n as listed in table [5] 
The species dependence is encapsulated in the 
function 



1 



3 = 1 



C (Z h x r ) + C 1 (Z i ,x I )T ] Zj>l, 



where 



C (^j,x r ) 
Ci(Zj,x r ) 



a (a+ 
&o(a;r) 



- ai(a; r ) lnZ^ 
6i(x r 



ln^-. 



We fit the coefficients ao(x T ), ai(x r ), 6o(a; r ), ^i(^r) 
piecewise over disjoint and x r ranges as follows: 



For 10 < Tj < 30 



0.973587- 0.0790899 In 



x x 



< 10 



a (x T ) 



ai(x T ) 

b (x r ) 
bi{x r ) 



0.806877 
0.652748 - 

0.766772 
-0.029204 + 
0.0853907 



x t > 10 
0.117973 \nx r 



Xr < 2.5 



x r > 2.5 

0.0367777 exp(-0.933151nx r ) 
0.053054 exp(-0.885854 lnx r ), 



whereas for 30 < F,-: 



a (x T ) 



ai(x T ) 



b (x T ) = 



bl(Xr) 



0.979517- 0.1279 In x r 
+0.0252858(lna; r ) 2 - 0.00204823(lnav) 3 
1.02454 + 0.0577405 lna; r x r < 5 



5. Application of the MLMR to obtain 
thermal and electrical conductivities of 
superburst progenitors 

Figure (fl"3)) shows the thermal conductivity cal- 
culated from numerical MD simulations in this pa- 
per and accurately reproduced by using the mix- 
ing rule prescription of Eq.(|40|). The conductivity 
of the amorphous outer crust is lower by an or- 
(46jfier of magnitude from estimates obtained using 
the impurity parameter formalism, Eq. (|2"5"|) . Fig- 
ure (JI3]) (see also Fig.©) also show that (Z 2 ), 
and as a consequence the conductivity, are less 
sensitive than Qimp to reaction processes that re- 
duce crust impurity abruptly in the outer crust, 
such as (1) electron captures on very proton-rich 
nuclei produced in the rp-process (and therefore 
far away from beta-stability at the top of the 
crust), (2) charge particle capture and fusion re- 
actions destroying lighter nuclei, and (3) (7,n) 
and (n, 7)-driven rearrangement of abundances in 
mass chains. 

We note that in the limit of a classical amor- 
phous solid, the Coulomb logarithms for elec- 
trical and thermal conductivity are equal, i.e. 
lnA CT = lnA K , therefore a and k are related by 
the Wiedemann- Franz law, Eq. (fT2|) . Thus while 
we only show a plot of the outer crust thermal 
conductivity as function of density in fig. (fl"3|) . the 
composition-dependent electrical conductivity is 
related very simply to it. 

X-ray burst ash composit ion differences will ex- 
ist between 1-zone models (jSchatz et al.l 1200 II ) of 
an X-ray Burst as op posed to a multi-zone model 



1.09301 + 0.0135052 lna; r x 
0.00841864- 0.0124097 In x r 

-0.0193203 x T > 10 
0.0225627 + 0.0242463 lnx r 

0.0422014 + 0.0117571 lnx r 

-0.0193203 10 < x T 



> 5 



(jWooslev et al.ll2004f) . Our purpose here is not 



to study all possible compositions in a superburst 
progenitor, but rather to show the applicability of 
the Microscopic Linear Mixing Rule to an arbi- 
trarily complex composition with the heterogene- 
ity expected of X-ray burst ashes. Once compo- 
sition is known, MD permits us to determine the 
5 < x r < 10 physical state and the conductivity of crust mat- 
ter. Since the thermal conductivity is a critical 



< 10 



< 5 
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ingredient setting the crust thermal gradient, an 
uncertainty in this quantity makes it difficult to 
predict whether shallow crust conditions ~ 10 9 
g/cc are conducive to superburst ignition. We 
will leave the full thermal profile calculation to 
a later work, but we point out here that a conduc- 
tivity significantly lowered by impurities, such as 
the MD results show for the outer crust which is 
in the amourphous outer crust can foster a steep 
thermal gradient suc h as would be required for su- 



perburst ignition. In iBrownl (|2004f ) a local stabil- 
ity analysis for accretion rate M — 3 ■ 10 17 g s _1 
shows that superburst ignition at column depth 
y w 10 12 g cm~2 is favored by an amorphous crust. 
Further, such a crust renders the thermal profile 
relatively insensitive to the core neutrino luminos- 
ity, whether from modified URCA processes (in- 
efficient core cooling) or enhanced due to direct 
URCA or pionic reactions (efficient core cooling). 
These results were confirmed by the glob al stabil- 
ity analysis o f lCooper fc Naravanl (|2005l) . 



A complete analysis of thermal structure re- 
quires an understanding of conductivity in the 
inner crust and therefore of nuclear processes at 
neutron drip which may reduce impurity signif- 
icantly (but not completely) from that of X-ray 
burst ashes. Wc arc currently performing these de- 
tailed studies and will publish the results shortly. 

We would like to thank Sanjay Rcddy and Dany 
Page for stimulating discussions. This work was 
performed for the U.S. Department of Energy by 
Los Alamos National Laboratory under contract 
DE-AC52-06NA25396. 
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A. Details of MD simulations 

Some technical details of the MD simulations are collected in table [3j namely 

• N sp = number of nuclear species; 

• N = total number of particles in the MD simulations; 

• N e q = number of time steps used initially to let the system equilibrate; 

• -^Vrun = number time steps of the MD run after equilibration; 

• = total number of wave-vector norms used to calculated S zz (k) with a\k < 22, see Eqs. (|34j|3"5")) . 

In all the simulations, the time step St is chosen such that St = 0.01/o-> p , where lu p is the mean ion plasma 
frequency. The total energy is very well conserved during the simulation (better than one part in 10 7 ). 



B. Relativistic dielectric function 



The Random-Phase A pproximation static dielectric function of the relativistic electron gas in its ground 
state is ( Jancovici Il962h 



+ x 

x 2 smh x r + \J 1 + x 2 — — In 



6x 2 x 



1 -x 



1 - 2x^x 2 
6x 2 x 



yl + x 2 x 2 In 



y/1 + X 2 X 2 +X^1 +X 2 , 



y/T+ x 2 x 2 — xy 1 + x 2 



(Bl) 
(B2) 



where x = fc/2fcp and /ctf is the (nonrclativistic) Thomas- Fermi (TF) inverse screening length 



nT 



ey/12nm e n e 
Hkp 



(B3) 



C. Additional structure factors 

Fig. (|14[) and (jT5]) shows the charge-charge structure factors S zz (k) and their mixing rule approximation 
for some of the compositions of table [2] but temperatures different than 0.5 GK, as listed in table [31 



17 



REFERENCES 

Baiko, D.A., Kaminker, A.D., Potekhin, A.Y., 
Yakovlev, G.G. 1998, Phys. Rev. Lett., 81, 5556 

Brown, E.F. 2000, ApJ, 531, 988 

Brown, E.F. 2004, ApJ, 614, L57 

Cackett, E. M. et al. 2006, MNRAS, 372, 479 

Cassisi, S., Potekhin, A.Y., Pictrinfcrni, A., Cate- 
lan., M., Salaris, M. 2007, ApJ, 661, 1094 

Chabrier, G., Ashcroft, N.W. 1990, Phys. Rev. A, 
42, 3384 

Cooper, R., Narayan, R. 2005, ApJ, 629, 422 

Cumming, A., Bildsten, L. 2001, ApJ, 559, L127 

Cusak, N., Endcrby, J.E. 1960, Proc. Phys. Soc., 
75, 395 

Daligault, J. 2006, Phys. Rev. Lett., 96, 065003 

Daligault, J., 2006, Phys. Rev. E, 73, 056407 

DeWitt, H.E., Slattery, W. 2003, Contrib. Plasma 
Phys., 43, 479 

Flowers, E., Itoh, N., 1976, ApJ, 206, 218 

Gupta, S., Brown, E.F., Schatz, H., Moller, P., 
Kratz, K.-L. 2007, ApJ, 662, 1188 

Gupta, S., Kawano, T., Moller, P. 2008, 
Phys. Rev. Lett., 101, 231101 

Haensel, P., Zdunik, J. L. 1990, A&A, 227, 431 

Haensel, P.,Kaminker, A. D., Yakovlev, D.G 1996, 
A&A, 314, 328 

Haensel, P., Zdunik, J. L. 2003, A&A, 404, L33 

Haensel, P., Zdunik, J. L. 2008, A&A, 480, 459 

Hamaguchi, S., Farouki, R.T., Dubin, D.H.E. 
1997, Phys. Rev. E, 56, 4671 

Horowitz, C.J., Caballero, O.L., Berry, D.K. 2009, 
Phys. Rev. E, 79, 026103 

Itoh, N., Hayashi, H., Kohyama, Y., 1993, ApJ, 
418, 405; 436, 418(E)(1994) 

Itoh, N., Kohyama, Y., 1993,ApJ, 404, 268 



Itoh, N., Hayashi, H. , Nishikawa, A., Kohyama, 
Y. 1996, ApJS, 102, 411 

Itoh, N., Uchida, S., Sakamoto, Y., Kohyama, Y., 
2008, ApJ, 677, 495 

Jancovici, B. 1962, Nuovo Cim., 25, 428 

Jones, P.B., 2005, Phys. Rev. D, 72, 083006 

Kaminker, A.D., Pethick, C.J., Potekhin, A.Y. , 
Thorsson, V., Yakovlev, D.G. 1999, A&A, 343, 
1009 

Kuulkcrs, E. ct al. 2004, AIP Conf. Proc., 714, 257 

Kuulkcrs, E. 2004, Nuc. Phys. B Proc. Supp., 132, 
466 

Mott, N.F., 1934, Proc. Roy. Soc. A, 146, 465 

Potekhin, AY., Baiko, D.A., Haensel, P., 
Yakovlev, D.G. 1999, A&A, 346, 345 

Potekhin, A.Y., Chabrier, G., Rogers, F.J. 2009, 
larXiv:0812.4344fr 2 

Rosenfcld, Y., 1993, Phys. Rev. E, 47, 2676 

Schatz, H., et al. 2001, Phys. Rev. Lett., 86, 3471 

Shternin, P.S. et al. 2007, MNRAS, 382, L43 

Shternin, P.S. and Yakovlev, D.G., 2009, 
Phys. Rev. D, 79, 123004 

Trunbull, D., Cecil, R.E., 1950, J. Appl. Phys., 21, 
804 

Young, DA., Corey, E.M., Dewitt, H.E. 1991, 
Phys. Rev. A, 44, 6508 

Wallace, D., 1997, Phys. Rev. E, 56, 1981 

Wijnands, R. 2001, ApJ, 554, L59 

Woosley, S. et al. 2004, ApJS151, 75 

Yakovlev D.G., 1987, Sov. Astron., 31, 347 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 



18 




Fig. 1. — Evolution with depth of dimcnsionless 
paramete rs obtained u s ing th e crustal composition 



profile of iGupta et alJ (|2007l ). r o ff is the effective 
ion coupling parameter defined in Eq.Q, a/as is 
the ratio of the mean interionic distance to the 
electronic Bohr radius, x T is the relativity param- 
eter, G = T/T-p is the degeneracy parameter. The 
vertical dashes indicate the compositions studied 
with MD and listed in table [2 The dots show the 
ratio between the impurity parameter Qi mp de- 
fined in Eq. ([M)l and (Z 2 ) . As discussed in section 
[5l (Z 2 ) is less sensitive than Q- imp to composition 
changes at densities ~ 10 7 , 2 x 10 s , 10 10 and 10 11 
g.cm -3 . 
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Fig. 2. — The impact parameter formalism illus- 
trated. (Top panel) A perfect crystal with a low 
concentration of impurities substituted at lattice 
sites can be approximated as being made of two 
uncorrelated parts: the perfect lattice together 
with the residual charge. The total electron- 
ion scattering cross-section is the sum of the 
Bloch electron-lattice scattering cross-section plus 
the sum of uncorrelated binary electron-residual 
impurity collsional cross-section. (Lower panel) 
As the concentration of impurities increases, this 
splitting is less and less justifiable. 
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Fig. 3. — (color) The structure factor S zz {k) and 
the corresponding integrand for the Coulomb log- 
arithm Eq. (TT4"|) with (red line) and without (green 
line) the relativistic correction. The upper panel 
corresponds to the composition #1 of table [5] at 
T = 0, 27 GK (r cfr = 60 and (Z) = 37). The lower 
panel was obtained with compostion #9 of table 
[lat T = 0.8 GK (r off = 153 and (Z) = 34). The 
vertical dashed lines mark the limit of integration 
2/cp in the Coulomb logarithm (fT4")) for liquid met- 
als with effective ionic charges Z c h = 1,2,3 and 
4. Its position relative to the main peak in S zz (k) 
clearly depends on the valency of the metal ion. 
In the crust of a neutron star, (Z) >> Z e e, and 
consequently the range of integration encompasses 
all first peaks in S zz (k). Physically, in neutron 
star crusts, large momentum transfers (close en- 
counters) significantly contribute to the Coulomb 
logarithm. 



Fig. 4. — (color) Contribution of the relativistic 
correction to the Coulomb logarithm (fT4")) in the 
region 10 7 < p < 10 10 g.cm -3 . lnA c ; corresponds 
to the nonrelativistic limit of Eq. (TT4"]) obtained by 

setting the bracketed term 1 — jj^i to 1- 
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100 




Fig. 5. — Composition profile at two crustal 
depths corresponding to mixture #1 (upper panel) 
and #12 (lower panel) of tablc[2] The figure shows 
the number concentration as a percentage as a 
function of the nuclear charge Z. At a given Z, the 
horizontal tics indicate different nuclear masses A. 
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Fig. 6. — (color) The structure factor S zz (k) ob- 
tained with MD simulations (red lines) and with 
the mixing rule formula (blue lines) for the 
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Fig. 9. — (color) (upper panel) Pair distribution 
function g{r) between particles of most abundant 
species of compositions #1, #3, #6 and #15 de- 
fined in table O (lower panel) Pair distribution 
function between the most abundant and another 
less abundant species for the same compositions of 
the upper panel. 
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Fig. 10. — (color) Cumulated integral InA(K) 



1 - 



l+x? \2k F J 

tained using the S zz (k) from 



J K dkk 3 



v(k) 
e(k) 



ob- 



MD simulations 
(full blues line) and from the mixing rule Eq. (l39)) 
(red dots). The Coulomb logarithm Eq. (fH|) is 
In A(2fcp). The upper and lower curves correspond 
to the mixture (1) and (9) defined in table (J5J), 
respectively. 
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Fig. 11. — OCP Coulomb logarithm as a function 
of the Coulomb coupling T for several values of the 
relativistic parameter x T . 
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Fig. 12.— OCP Coulomb logarithm lnA OCP at 
x T = 1 and x T = 10. The red solid lines arc 
obtained from Eq.(??) with the OCP structure 
factors of lYoung et al.l (|1991l ) valid for 5 < T < 
225. The blue dotted lines are from our pre- 
scription Eq .d42l) and fit the values obtained with 



(IH £ 

(1991 



Young et alJ (|l99lh very well and smoothly extend 



into the T > 225 regime. 
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Fig. 13. — Thermal conductivity as a function of 
density in the crust of an accreting neutron star 
with M w O.lAfEdd- The molecular dynamics 
(MD) calculations (squares) which are accurately 
reproduced by the linear mixing rule (j40|) (red line) 
are contrasted with the electron-impurity scatter- 
ing conductivity as calculated using Eq. (|25|) over a 
range of four orders of magnitude in crust density. 
As discussed in sectional (Z 2 ), and as a conse- 
quence the conductivity, is less sensitive than Q lmv 
to composition changes at densities ~ 10 7 , 2 x 10 8 , 
10 10 and 10 11 g.cm" 3 . 
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Fig. 14. — (color) The structure factor S zz (k) ob- 
tained with MD simulations (blue lines) and with 
the mixing rule formula (|39|) (red line) for the mix- 
tures #1-4 defined in table SI 
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Fig. 15. — (color) Same as figure [T4l for composi- 
tion #5-8 listed in table H 
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Table 1 

Entropy variation As = AS/nVksT on melting (i.e. from BCC to fluid) of the OCP as a 

FUNCTION OF THE INVERSE DIMENSIONLESS SCREENING LENGTH K = ak sc . 



K 


r 


As 





171.8 


0.8516 


0.2 


173.5 


0.8526 


0.4 


178.6 


0.8551 


0.6 


187.1 


0.8580 


0.8 


199.6 


0.8599 


1. 


217.4 


0.8591 



Note. — r m is the 
coupling paramater 
Tocp at melting. 
The data were cal- 
culated using the 
equa tion of state 



of 


Hamaguchi et al. 


(1997 


)■ 



Table 2 

List of the mixtures considered in this work extracted from Igupta et alI ( 20071 ). 



mixture 


sp 


N sp 


P (g-cm 3 ) 


E F (MeV) 


akp 


r off = (z 5 / 3 )r e 


(Z) 


(Z 2 ) 


1 


450 


55 


8.59 10 6 


0.45 


6.44 


32.3 


37 


1597 


2 


419 


50 


1.01 10 8 


1.38 


6.46 


71.3 


37 


1537 


3 


451 


52 


2.49 10 8 


2.00 


6.44 


97.9 


38 


1564 


4 


432 


52 


4.33 10 8 


2.49 


6.52 


118.4 


38 


1572 


5 


405 


50 


7.08 10 8 


3.00 


6.42 


137.4 


38 


1549 


6 


400 


53 


1.00 10 9 


3.42 


6.41 


153.2 


38 


1536 


7 


398 


53 


2.87 10 9 


5.00 


6.37 


179.5 


35 


1359 


8 


406 


52 


4.85 10 9 


6.01 


6.31 


236.4 


36 


1401 


9 


411 


49 


1.00 10 10 


7.72 


6.26 


263.1 


34 


1275 


10 


404 


49 


2.17 10 10 


10.0 


6.20 


343.8 


34 


1171 


11 


372 


52 


3.74 10 10 


12.0 


6.18 


390.6 


33 


1113 


12 


355 


50 


4.26 10 10 


12.5 


6.18 


399.9 


33 


1090 


13 


363 


54 


4.81 10 10 


13.0 


6.14 


410.6 


32 


1073 


14 


341 


50 


5.99 10 10 


14.0 


6.15 


436.8 


32 


1119 


15 


344 


49 


7.37 10 10 


15.0 


6.15 


461.7 


32 


1103 


16 


341 


46 


9.13 10 10 


16.0 


6.04 


468.71 


31 


1041 


17 


206 


24 


1.75 10 11 


19.7 


5.95 


538.8 


30 


964 



Note. — N°* 1S is the original number of species in iGupta et all (|2007h . N sp is the number of 
species used in the MD simulations, p is the mass density, Ep the Fermi energy, T the temperature, 
akTF the dimensionlcss TF screening length where a is the interparticle distance given by Eq.(??), 
ah? is the dimensionlcss Fermi wavevector, T is the coupling parameter defined by Eq.Q. 



Table 3 

Simulation parameters used for the mixtures defined in table El 



mixture 


N sp 


N 




iVrun 


N k 


i 
i 


•) 


071 & 


1 

1U 


1 

1U 


1 AA^\ 


n 
z 


u 


1 QC«/L 


1 

1U 


1 

1U 


Z^u^i 


9 


•JZ 


lyooi 


1 n 5 

1U 


1 

1U 


z^y * 


4 






1U 


1U 


9^0*3 
zOUo 


r. 



u 


lyooi 


1U 


1U 


9/1Q7 
Z4y f 




o o 


14oOZ 


1 n 5 

1U 


1 n 5 

1U 


91 07 
Z±U ( 


7 


o o 


1 QC/IQ 


1 n 5 

1U 


1 n 5 

1U 


9^17 
ZOl / 


o 


z 


1 QQK7 


1 n 5 

1U 


1 n 5 

1U 


9^17 


9 


49 


14847 


10 5 


10 5 


2107 


10 


49 


19859 


10 5 


10 5 


2506 


11 


52 


19870 


10 5 


10 5 


2869 


12 


50 


19875 


10 5 


10 5 


3095 


13 


54 


19871 


10 5 


10 5 


3287 


14 


50 


19898 


10 5 


10 5 


3350 


15 


49 


29888 


10 5 


10 5 


4299 


16 


40 


39903 


10 5 


10 5 


4972 


17 


24 


19937 


10 5 


10 5 


2143 



Table 4 

Compositions used for the MD calculations shown in FigsIT41 and [T5l 



mixture 


T (GK) 


r cff = (z 5 / 3 )r e 


1 


0.27 


GO 


2 


0.73 


49 


3 


0.81 


GO 


4 


0.84 


70 


■5 


0.85 


81 


6 


0.85 


90 


7 


0.84 


106 


8 


0.83 


142 



Note. — For a given mixture, the com- 
position is the same as the one given in 
table [TJ The only different is the temper- 
ature, and thcrfore T e g. 



Table 5 

Fitting parameters a n and b n , n = 1,2, 3, to be used in Eqs.(|44|) and l4"5l 



n a n b n 



1 48.12 1.95 

2 29.23 11.74 

3 0.77 92.18 



